The Exponential Time Complexity of Computing 
the Probability That a Graph is Connected* 

Thore Husfeldt * Nina Taslaman * 

Abstract 

£SJ ' We show that for every probability p with < p < 1, computation 

of all-terminal graph reliability with edge failure probability p requires 
O ■ time exponential in 0(m/ log 2 m) for simple graphs of m edges under the 

' Exponential Time Hypothesis. 

1 Introduction 

• Graph reliability is a simple mathematical model of connectedness in networks 

| that are subject to random failure of its communication channels. This type 

of stochastic networks arise naturally in, e.g., communication or traffic control; 
| see pQ for an extensive survey of application areas. 

For a connected graph G = (V, E) and probability p, the all-terminal relia- 
bility R(G;p) is the probability that there is a path of operational edges between 
every pair of nodes, given that every edge of the graph fails independently with 
I probability p. For example, with p = g, the all-terminal reliability of the graph 

^sO . <£j? is , and the all-terminal reliability of the graph hC is 

m 

O; R(G; P )= V pl^l(l-p)W. (1) 
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In general, for a connected, undirected graph G = (V, E) with edge-failure 
probability p, the all-terminal reliability can be given as 



spanning 
connected 



For example §) = (|) 5 + 5 • § ■ (f ) 4 = ±±§. Computing i?(G;p) directly 

from ([T]) can take up to 0(2 m ) operations, where m is the number of edges 
in G. An algorithm of Buzacott |1] solves the problem in time 3 n n°^ for 
graphs of n vertices. Further improvements exist [2J, but remain exponential in 
n; this is explained by the exp(Q(n.)) lower bound of [5]. On the other hand, 
subexponential time algorithms have been found for some restricted classes of 
graphs. For example, the problem can be solved in time exp(0(y / n)) for planar 
graphs 14 . A natural question is then whether the complexity can be reduced 
for for further classes of graphs. Especially, from an applications point of view, 
the case of simple graphs is interesting. 
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It is clear that R(G; 0) = 1 and R(G; 1) = for all connected graphs, so 
for some values of p the problem is trivial. One may ask what the situation is 
for values close to these extremes. Moreover, for p = I it is easily seen that 
(P) equals the number of connected, spanning subgraphs of G, divided by 2™. 
This is an interesting enumeration problem in itself, and one could be tempted 
to hope for a better algorithm than exp(0(n)), because a related enumeration 
problem, the number of spanning trees of G, can be solved in polynomial time 
by Kirchhoff's matrix-tree theorem 

Result 

We give a lower bound on the problem of computing all-terminal graph reliabil- 
ity for the class of simple graphs for all nontrivial p, in the framework recently 
proposed by Dell et al. [5]. In particular, we work under the counting exponen- 
tial time hypothesis: 

(#eth) There is a constant c > such that no deterministic algorithm can 
compute ^3-Sat in time exp(cn). 

This is a relaxation of the exponential time hypothesis (eth) of Impagliazzo 
et al. 9 , so our results hold under eth as well. The best current bound for 
#3-Sat is 0(1.6423") Q3]. 

Theorem 1. For any fixed probability p with < p < 1, computing the all- 
terminal reliability R(G;p) of a given simple graph G of m edges requires time 
exponential in f2(m/ log 2 m) under #ETH. 

In particular, the bound holds for p = h, i.e., counting the number of con- 
nected spanning subgraphs of a given graph. 

We have expressed the lower bound in terms of the parameter m, the number 
of edges of the input graph. Since n < m for connected graphs, the result implies 
the lower bound exp(f2(ro/ log 2 n)) in terms of the parameter n, the number of 
vertices of the input graph. Moreover, the fi(m/ log 2 m) lower bound together 
with the exp(0(n)) algorithm from [H [2] shows that the hard instances have 
roughly linear density, ruling out a better algorithm than exp(0(n/ log 2 n)) also 
for the restricted case of sparse graphs. 

Our bound does not quite match the best known upper bound exp(0(n)) of 
[4l [2] . This situation is similar to the bounds reported in [5] for related problems 
on simple graphs, which also fall a few logarithmic factors (in the exponent) 
short of the best known algorithms. The bound does, however, suffice to separate 
the complexity of reliability computation from the exp(0(y / n)) bound for the 
planar case [T4] . 

Graph polynomials 

Expression (fTJ, viewed as a function p M> R(G;p) for fixed G, is known as 
the reliability polynomial of G, an object studied in algebraic graph theory [6j 
Sec. 15.8]. For example, R{<&;p) = (1 - pf + bp(l~p) A . 

Arguably, the most important graph polynomial is the bivariate Tutte poly- 
nomial T(G;x,y), which encodes numerous combinatorial parameters of the 
input graph G, and whose restriction to certain lines and curves in the xy-plane 
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Hexp(fi(n)) [5] 

I | exp(f2(m/ log 2 m)) (this paper) 

□ exp(n(m/log 3 m)) [5] 

□ unless P = #P [TO] 
■ n°W HQ] 



Figure 1: Exponential time complexity under #eth of the Tutte plane for simple 
graphs. 

specialize to other well known graph polynomials. The reliability polynomial is 
essentially a restriction of this polynomial to the ray { (1,2/) : y > 1 }. The com- 
plexity of computing T(G; a:, y) at various points (x, y) with respect to an input 
graph G is very well-studied in various models of computation, and the present 
paper thus establishes lower bounds for simple graphs along the mentioned ray, 
which was left open in a recent study [5] to completely map the exponential 
time complexity of the "Tutte plane" . 

Related work 

The structural complexity of all-terminal graph reliability was studied by Provan 
and Ball jl3j . For any probability p, with < p < 1, it is shown that computing 
R(G;p) for given G is hard for Valiant's counting class #P [TBJ- The reductions 
in [T3] do not preserve the parameters n and m, so that the running time bounds 
under #ETII implicitly provided by their techniques are typically exponential in 
Vl{n 1 l k ) for some fc. 

The reliability problem under consideration in this paper admits a number 
of natural extensions. 

1. We can consider the computational problem of finding the reliability poly- 
nomial itself, instead of its value at a fixed point p. The input to this 
problem is a graph, and the output is a list of coefficients. For example, 
on input <£p the output should give i?(<£j? ;p) = 4p a — 15p 4 + 20p 3 — 10p 2 + l. 

2. We can associate individual probabilities to every edge. For example, the 
graph \ \ becomes disconnected with probability |. 

3. We can consider multigraphs like ° — =o, but with the same edge weight p. 
As indicated by the examples (for p = i), the multigraph case is a special 
case of the individually edge- weighted case, a fact that we will use later. 

All of these problems are at least as hard as the problem under consideration 
in the present paper. Lower bounds of size exp(f2(m)) are given in [S] or follow 
relatively easily (see H2. 3p . 

A recent paper of Hoffman 8 studies the complexity of another graph poly- 
nomial, the independent set polynomial, in the same framework. 
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2 Preliminaries 



We will only be concerned with undirected graphs. For a graph G — (V, E) let n 
denote the number of vertices, m the number of edges, and for any subset A C E 
let k(A) denote the number of connected components in the subgraph (V, A) 
(especially, k(A) = 1 means that the edge subset A is spanning and connected). 
Also, for graph polynomials P and Q, we write P(G;x) ~ Q(G';x') if the two 
expressions are equal up to an easily computable factor. 

2.1 Weighted reliability 

The reliability polynomial can be formulated as a restriction of the Tutte poly- 
nomial, which for an undirected graph G is given by 

T(G; x,y)=J2( x - ir {A) - K(E) (y - 1)< A )+W~\ V \ . 

ACE 

Note that k(E) = 1 in our case. We find the reliability polynomial along the 
ray { (1, y) : y > 1 } in the so called Tutte plane, as R{G;p) ~ T(G; 1, in 
full detail: 

i?(G;p)=p m -"+ 1 (l-p)"- 1 T(G;l,l/p), (0<p<l). (2) 

For complexity analysis of the Tutte polynomial, it has proved a considerable 
technical simplification to consider Sokal's multivariate Tutte polynomial |15) . 
Here the graph is equipped with some weight function w: E — > K, and the 
polynomial is given by 

Z{G;q,w)= £ w(A)q"W , 

ACE 

where w(A) — YiecA w ( e ) ^ s * ne edge-weight product of the subset A. For 
constant edge weights w — y — 1 we have Z(G; q, w) ~ T(G; x, y) with g = 
(x — l)(y — 1). The "reliability line" x = 1 in the Tutte plane thus corresponds 
to q = in the weighted setting, where Z vanishes, so instead we will consider 
the slightly modified polynomial 

Z(G;q,w)=q- 1 Z(G;q,w). 

At q — 0, this gives a weighted version of the reliability polynomial: 

Definition 1. (Weighted reliability polynomial) 

For a connected, undirected graph G = (V, £"), the weighted reliability polyno- 
mial of G is given by 

R(G;w) = Z(G;q,w)\ q=0 = £ w(A) . (3) 

k(A) = 1 

For constant edge weight w > we have -R(G; w) = w n ~ 1 T{G; 1,1 + w), so for 
< p < 1 we can recover the reliability polynomial through as 

R(G;p)=p m R(G;l/ P -l). (4) 
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2.2 Graph transformations 



A classical technique for investigating the complexity of the Tutte polynomial at 
a certain point (a/, y') of the Tutte plane, is to relate it to some already settled 
point (a;, y) via a graph transformation <p, such that T(G; x', y') ~ T(ip(G); x, y). 
For the weighted setting we have the following rules, which are simple general- 
izations of I7J Sec. 4.3]. (See Appendix lA.il ) For a graph G = (V,E) with edge 
weights given by w: 

Lemma 1. If ip(G) is obtained from G by replacing a single edge e G E with a 
simple path of k edges P — {e\, e^} with w(ej) = W{, then 

R(ip(G);w) = C P ■ R(G; w[e i— > w'}) , 

where 

h 

111 1 . 

— - = 1 1 and Cp = — TT Wi . 

W W\ Wk w ± - L 



Lemma 2. If ip(G) is obtained from G by replacing a single edge e^E with a 
bundle of parallel edges B = {e\, . . . , e^} with w(ej) = Wi, then 

R(tp(G);w) =B(G;w[e4 w'}) , 

where 

k 

w' = -l + Y[(l + w t ). 



Corollary 1. If ip(G) is obtained from G by replacing a single edge eeE with 
a simple path of k edges of constant weight w, then 

R((p(G);w) = kw*- 1 ■ R(G; w[e ^ w/k]) , (5) 

and if it is obtained from G by replacing e € E with a bundle of k parallel edges 
of constant weight w, then 

R(<p(G);w) = i?(G;w[e h> (1 + w) k - 1]) . (6) 

These rules are transitive [3 Lem. 1], and so can be freely combined for 
more intricate weight shifts. To preserve constant weight functions we need to 
perform the same transformation to every edge of the graph. This calls for the 
graph theoretic version of Brylawski's tensor product for matroids [3]. We found 
the following terminology more intuitive for our setting: 

Definition 2. (Graph inflation) 

Let H be a 2-terminal undirected graph. For any undirected graph G — (V, E), 
an H -inflation of G, denoted G ® H , is obtained by replacing every edge xy € E 
by (a fresh copy of) H , identifying x with one of the terminals of H and y with 
the other0 

If H is a simple path of k edges, G ® H gives the k-stretch of G. Similarly, a 
bundle of k parallel edges results in a k-thickening, of G. 

1 This can, in general, be done in two different ways, resulting in graphs that need not be 
isomorphic. However, the Tutte polynomial is blind to this difference. See extensive footnote 
in [5], Section 5.1. 
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2.3 Hardness of computing coefficients 

Our pivot for proving Theorem[T]will be the following hardness result, which says 
that even when restricted to fixed hyperbolas (x — l)(y — 1) = q, computing the 
full Tutte polynomial is hard. This is an extension to the case q = of Lemma 2 
in [S], and the proof is given in Appendix 2. 

Lemma 3. Under #eth, computing the coefficients of the polynomial w i— > 
Z(G;q,w) for given simple graph G and rational number q ^ {1,2} requires 
time exponential in f2(m). 

Since R(G;p) is essentially R(G;w) restricted to positive constant weight 
functions, and since R(G; w) = Z(G;0,w), the following is immediate: 

Corollary 2. Under #eth, it requires time exponential in f2(m) to compute 
the coefficients of the reliability polynomial. 

3 Bounce graphs 

As a first step towards Theorem [TJ we present here a class of graph transfor- 
mations whose corresponding weight shifts for the reliability polynomial are 
all distinct. These transformations are mildly inspired by fc-byte numbers, in 
the sense that each has associated to it a sequence of length k, such that the 
lexicographic order of these sequences determines the numerical order of the 
corresponding (shifted) weights. Each transformation is a bounce inflation: 

Definition 3. (Bounce graph) 

For positive numbers h (height) and I (length), the {h,l)-bounce is the graph 
obtained by identifying all the left and all the right endpoints of h simple paths 
of length I. Given a bounce sequence, S — (s\, S2, . . . , s/c), of k numbers Si > 1, 
the corresponding bounce graph, B$, is the (simple) graph obtained by concate- 
nating k (h, Z)-bounces by their endpoints, where the height starts at 1 for the 
first bounce and then increases by one for each follower, and the length of the 
ith bounce is Si. 



The length of a bounce graph is the number of bounces in it (or, equivalently, 
the height of the highest bounce). 

Inflation by a bounce graph has the following weight-shifting effect, from i?'s 
perspective: 

Lemma 4. For a graph G, bounce sequence S — (si, S2, . . . , s/J and w > 0: 




(4, 2)-bounce 



S = (3, 2,3,2) 



G 



where 

k k 

— = E n I /\i i and C S = — •n^- 1 ) i (( w + Si ) i - s i) . 
»s r-f (1 + w/sj) 1 - 1 «"s 1J ; 

Z=l 2=1 



Proof. Starting out with G ® -Bs, we will look at the effect of replacing one 
of the m copies of Bs with a single edge e. We show that, with ip denoting this 
operation: 

R{G ® S s ; w) = C s - %(G ® w[e ^ , (7) 

where has the above form, and w has the old value w on all unaffected edges. 
The lemma follows from performing ip for every copy of Bs in G® B$- 

The first step towards transforming a bounce graph (say, o-^^*V^\) into a 
single edge, is to replace each path of each bounce in it by a single edge. Apply- 
ing J5]) of Corollary Q] to each path of the ith bounce gives a factor (si-u;^ 1 ) 1 to 
the polynomial, and each edge in the resulting (h, l)-bounce gets weight w/si in 
the modified graph. Repeating this process for every bounce gives a simplified 
bounce graph (»— ~^^\) in a transformed graph </>(G ® Sg) such that 

i?(G®B 5 ;^)= fllC^" 1 )^ •%(G®S s );w') ) 

with w' taking the value io/sj for every edge in the ith bounce of the simplified 
bounce graph, and the old value w outside it. Next, we successively replace 
each of its (h, l)-bounces by a single edge, to get a simple path („ — o — <• — o) of 
length k (with non-constant edge weights). From (|6]) of Corollary [TJ we know 
that this does not produce any new factors for the polynomial, but the weight 
of the ith edge in this path will be given by 

Wi = (1 + w/si) 1 - 1 . 

Finally, we compress the path into a single edge e. A single application of 
Lemma [1] then gives the result in (JJJ. □ 

If Lemma [3] is our pivot for proving Theorem [TJ the following result is the 
lever: 

Lemma 5. For any size m, there exist m + 1 distinct, simple bounce graphs Bs 
of size 0(log 2 m), such that for any two associated bounce sequences S and T : 

S >i ex T =$> w s <w T (8) 

for all w > 6. 

Proof. The set of bounce sequences S = (s\, . . . , si) of length I — log(m+l) 
and with each Sj € {2,3}, provides m + 1 different simple bounce graphs of 
promised size. To show that any two of them satisfy equation we will look 
at the difference 

A / N 1 1 

As,t(w) = , 

Ws WT 
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and show that As,t(w) > for w > 6. 

Let k be the first index where the sequences differ, say Sk = 3 and tk = 2. 
We then have 

1 1 1 

AS ' T{W) = (1 + W3)*-1 + ? (I 



(l + w/3) fe -l (l + u-/ Sl ) l -l 



1 - 

>/2) fe - 1 ~ ^ (T 



(l + ™/2)fc-l ^(l + ^/tfcJi-l 



This would be minimal if = 2 and t j = 3 for all i > fc, i.e. if 

A s ,rH = /(I + ifl/3) - /(l + w/2) , 

where 

If we could show that /'(x) < for x > xo then it would follow (e.g. from the 
mean value theorem) that f(x) > f(y) for x < y above xq. In particular, with 
xq = 3 this would prove our claim. To see that this is indeed the case, we look 
at the derivative 

f( x ) = 5? + V — 

J { ' (x k -l) 2 (^-l) 2 ' 

A bit of manipulation shows that the terms of the sum, let us call them Tj, 
satisfy 7$ > 2T i+ i for x > 3, so 

for x > 3, and we are done by the above argument . □ 



4 Evaluating the reliability polynomial is hard 

We are ready to prove Theorem [TJ We introduce the following notation for the 
problem of evaluating a graph polynomial P(G; x) at a given point x 

P-val(x): 

input A simple, connected, undirected graph G 
output The value of P(G; x) (a rational number) 

In this notation, the computational problem described in Theorem [JJ can be 
written as R-VAh(p), and the corresponding problem for weighted reliability 
with constant edge weight w is R-val(w). We will prove Theorem[JJby reducing 
the problem of computing coefficients of the polynomial p i-> R(G;p), to the 
problem R-VAh(p) for any arbitrary fixed probability p with < p < 1. 
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Proof (of Theorem [ip. Let G be a simple graph with n vertices and m 
edges. We prove that R-val(w) requires time exponential in Q(m/ log 2 m) for 
any w > 0, which by Q gives the same bound for R-VAh(j>) for any p with 
<p < 1. 

Suppose we have an algorithm for R-val(w) for some fixed w > 0. For 
simplicity of exposition, first assume w > 6. From Lemma [3] we can easily 
construct m + 1 bounce graphs B$ such that each R(G Cg> -Bs; u;) gives us the 
value of R(G; ws) at some new weight ws, with all such ws distinct. Computing 
R-val(w) for each of these m + 1 bounce inflations, we get the value of R(G; w) 
at m + 1 distinct upvalues. Since the degree of this polynomial is m, that 
gives us the coefficients by interpolation. By Corollary [2 the whole process 
must then require time exp(51(m)). By Lemma [SJ each inflation G<g>B>s will 
have 0(m log 2 m) edges. Thus, for graphs of size 0{m log 2 m), the problem 
R-val(w) has a lower bound of exp(f2(m)). The claimed bound then follows 
from the fact that 

tp(m) = m log 2 m =>■ y> -1 (?7i) G O (mj log 2 m) . 

We turn to the case < w < 6. Given such a «;, choose a number k such 
that 

to' := (w/2 + l) k - 1 > 6. 

From the above, we know that for any simple graph G it takes time exponential 
in Q(m/ log 2 m) to evaluate R(G;w'). Now consider the fc-thickening-2-stretch 
of G, let us call this G'. This will be a simple graph of 0(km) edges, and 
from Corollary [1] it follows that the value of R(G';w) would give us the value 
of R(G; w'). Thus, computing the former must also require time exponential in 
Q(m/ log 2 m). Since G' has 0{km) edges, this gives a lower bound for _R-VAL(u;) 
exponential in VI (m/ (k log 2 m)) , which is J7(m/ log 2 m) as a function of m. □ 

5 Remarks 

For the multivariate Tutte polynomial Z(G\ q,w), the current lower bound for 
simple graphs is exp(0(m/ log 3 m)) 5j. One might ask whether the bounce 
graph construction could be used to improve this. The weight shift correspond- 
ing to a bounce inflation is in this case given by (for and q (£ {0, —2w}) 

-g- = TTf „ / / v l) ■ where r. i= fl + 

w S y \{(q/n) + iy-l ) 1 \ w) 

Unlike the expression in Lemma H] this is not a sum of powers, so the 'fe-byte 
number' analogy is remote, and there is now also a dependency on q which seems 
to make it difficult proving something like Lemma [5] under the same constraints. 
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A Supplementary proofs 



A. 1 Lemma Q] and [2] 

Proof (of Lemma[l] and\^). We repeat the arguments from the proof in [101 
Sec. 4. 3]: Let S be the set of subsets A C E \{e} that already span the whole 
graph G, i.e. 

S = {ACE\{e} : k{A) = 1} , 
and let T be the set of subsets that need the edge e to span the graph: 

T = {A C E \ {e} : k(A) = 2 and k(A U {e}) = 1} . 

With w' denoting the weight of the edge e in the original graph G, ([3]) gives, 
for both lemmas, 



R(G; w[e i— > to' 

]) = E w (^)(! + ">') + E W ( A W ■ 



(9) 



Ags 



We will compare the partial sums here to the corresponding ones obtained when 
we alter the graph. When ip is the operation described in Lemma [TJ we have 
(with P the set of edges in the path) 

i%(G) ;w) = J2 MA) ( w(P) + 5>(P\^))+E w(A)w(P) = 



Ags 



AGT 

fe 



/ fe fe \ 

x;w(A) n^+En^ +E w ( A )iT w - 

AGS \i=l 3=1 / AGT i=l 

Comparing corresponding sums to ([9]), it is easy to check that the expressions 
for to' and C p in Lemma [1] indeed make R(ip(G);w) = Cp ■ R(G;w[e H> to']). 

When ip is the operation described in Lemma [H we have (with B the set of 
edges in the bundle) 



£(¥>(G);w) = £w(A) 



AGS 



1+ £ w(A') 



V 



A'CB 
A'#0 



E w ( A ) 



/ 



AGT 



\ 



E w ( A ') 

,A'CB ; 
\AV0 / 



and Lemma [5] follows since 



]T w(A')=I]KH)-l. 



A'CB 
A'#0 



□ 
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A. 2 Deletion/contraction and Lemma [3] 



For an edge e = xy G E, let G \ e be the graph obtained by deleting e, and 
let G/e be the (multi)graph obtained by contracting e, i.e. by identifying the 
end vertices x and y before removing e . With w e denoting the weight of edge 
e, we have the following deletion/contraction reduction for the weighted Tutte 
polynomial (see [HI Sec. 4.3]) 

Z{G- q, w) = Z(G \ e; q, w) + w e ■ Z{G/e; q, w) . 

Note that if e G E is a bridge, then G \ e has one more component than G, 
while in any other case both G\e and G/e have the same number of connected 
components as G. Using the above identity and the fact that Z(G;q,~w) — 
q~ 1 Z(G;q,w), this gives: 

^^ G . qw s = \l-Z{G\e;q,w) + w e -Z(G/e;q,w) if e is a bridge, _ 
1 Z(G \ e; q, w) + w e ■ Z(G/e; q, w) otherwise. 



Proof ( of Lemma\Bj). The case q ^ is treated in Lemma 2 of [5] . For q = 0, 
we give a reduction to Lemma 1 from the same paper. Under our assumptions 
it says that w i-> Z(G; 0, w) cannot be computed faster than exp(fi(m)), where 
w is given by (for some set T of three edges): 

w(e) = <^ . (11) 

otherwise. 

The proof of this actually uses the restriction that G' — (V, E\T) is connected, 
so we can assume that this is the case. Thus, no edge in T is a bridge. Three 
applications of (fTU|) . to delete/contract these edges, gives 

Z(G;0,w)= Y, (-l) lCl Z(G c ;0,w), 

CC{1,2,3} 

for some graphs Gc of constant edge- weight w. These Gc's may contain loops 
and multiple edges from the contractions. To address this we look at the 2- 
stretch Gc <£> P3 of each Gc- This will give simple graphs of constant weight 
functions, and m applications of ([5]) from Corollary [T] gives 

Z(G C ® P 3 ; 0, w) = (2 W )™Z(G C ; 0, w/2) . 

If an algorithm could compute the coefficients of w h-> Z(G;0,w) faster than 
exp(f2(m)) for any simple graph G, it could be used to compute eight polyno- 
mials w 1— >• Z(Gc ® P3; 0, w) (one for each subset G). This would give us first 
w i-> Z(Gc]0, w) and then w h-> Z(G;0,w), faster than exp(fi(m)), which is 
impossible according to Lemma 1 from [5]. □ 
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